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ABSTRACT 

Patients  at  Boston's  Children's  Hospital  who  were  diagnosed 
as  having  Cerebral  Palsy  were  filmed  walking.   These  films  were 
digitized  and  translated  into  measurements  associated  with  leg 
motion.   Using  a  kth  nearest  neighbour  clustering  procedure,  five 
subpopulations  are  identified  from  the  modes  of  a  uniformly  consis- 
tent estimate  of  the  underlying  measurement-space  density. 

Graphical  profiles  and  F-ratios  were  used  to  identify  indi- 
vidual measurements  which  are  most  useful  in  distinguishing  the 
membership  of  the  various  observed  clusters .   Walking  velocity  and 
patterns  of  hip  and  ankle  movements  are  seen  to  be  markedly  differ- 
ent for  each  cluster.   Traits  associated  with  patients  in  two  of 
the  clusters  are  identified  with  severe  gait  pathologies,  while 
the  remaining  three  clusters  exhibit  characteristics  more  closely 
approximating  normal  gait . 

Using  the  recursive  partitioning  classification  rule,  it  is 
demonstrated  that  the  effects  of  corrective  measures  such  as  sur- 
gery or  assistive  devices  on  the  gait  patterns  of  individual 
patients  can  be  examined. 

KEYWORDS:  Cerebral  Palsy;  gait  analysis;  kth  nearest  neighbour 
clustering  procedures;  graphical  profiles;  recursive 
partitioning  classification  rule. 


1.   INTRODUCTION 
The  evaluation  of  human  movement  has  long  been  recognized  as 
a  necessary  and  important  aspect  of  the  clinical  workup  of  a  subject 
suffering  from  any  one  of  a  variety  of  disorders.   And  with  the 
numerous  developments  in  the  computer  and  electronic  technological 
applications  to  this  area  of  medicine  in  the  last  decade,  there  currently 
exist  a  plethora  of  instruments  and  methodologies  which  can  acquire, 
process,  and  display  human  movement  data  accurately  and  reproduceably 
within  a  clinically  useful,  short  period  of  time.   Although  such 
devices  and  methodologies  have  their  widest  clinical  application 
currently  in  evaluating  subjects  with  locomotive  disorder  (walking 
on  level  ground) ,  there  is  increasing  interest  and  every  reason  to 
believe  that  its  utilization  will  be  expanded  to  include  a  wide  vari- 
ety of  human  tasks;  for  example,  bending,  lifting,  going  up  and  down 
stairs  or  ramps,  writing,  feeding  oneself,  and  in  general,  performing 
all  the  activities  of  daily  as  well  as  recreational  living.   It  is 
also  clear  that  such  information  will  be  invaluable  to  a  wide  variety 
of  clinicians  working  in  the  fields  of  rehabilitation,  orthopedics, 
neurology,  and  sports  medicine,  to  name  but  a  few.   In  all  these  areas, 
testing  of  human  performance  could  be  utilized  to  (a)  understand  the 
pathodynamic  manifestations  of  a  given  organic  disorder,  and  (b)  monitor 
the  progress  of  an  untreated  as  well  as  a  treated  subject.   The  conse- 
quences that  this  could  have  towards  the  improvement  of  patient  care 
are  many. 
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I 
Unfortunately,  in  spite  of  the  continued  improvement  in  the  instru- 
mentation and  methodologies  of  evaluating  human  peformance,  objective 
measurements  of  human  performance  still  lack  the  widespread  everyday 
clinical  utility  that  X-ray  currently  has,  owing  to  the  inability  of 
the  clinician  to  properly  interpret  the  data.   The  major  reason  for 
this  is  that  the  evaluation  of  human  movement  requires  the  monitoring 
of  many  movements  and  many  forces  creating  these  movements,  each  one 
needing  to  be  monitored  many  times  in  a  given  time  period.   And  the 
problem  of  properly  interpreting  this  vast  amount  of  data  is  formid- 
able.  (See  however,  Sutherland  et  al  (1980),  for  a  statistical  study 
of  the  development  of  mature  gait.)   In  this  paper,  the  gait  measure- 
ments obtained  for  62  Cerebral  Palsy  patients  at  Boston's  Children's 
Hospital  are  used  to  illustrate  how  various  statistical  methodologies 
including  clustering  procedures  and  classification  rules  can  be  use- 
fully applied  to  the  analysis  and  interpretation  of  human  movement 
data. 

A  detailed  description  of  the  various  gait  measurements  obtained 
for  each  patient  in  this  study  is  given  in  Section  2.   In  Section  3,  a 
cluster  analysis  of  the  gait  measurements  data  is  presented.   The 
Cerebral  Palsy  patients  are  grouped  on  the  basis  of  the  various  gait 
measurements  using  the  kth  nearest  neighbour  clustering  procedure 
given  in  Wong  (1981) .   And  it  is  found  that  the  five  resulting  clusters 
can  be  identified  with  different  severeness  levels  of  abnormal  gait. 
Various  graphical  profiles  are  given  in  Section  4  to  illustrate  the  re- 
lative characteristics  of  these  clusters.   F-ratios  are  also  computed 
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and  are  used  to  identify  individual  variables  that  are  important  in  dis- 
tinguishing the  observed  clusters.   The  variables  identified  include 
walking  velocity,  cadence  of  steps,  and  step  length.   These  variables 
agree  with  the  measurements  cited  in  Sutherland  et  al  (1980)  as  impor- 
tant determinants  of  gait  development  in  normal  children,  and  they  are 
useful  in  illustrating  the  types  of  motions  which  are  most  affected  by 
Cerebral  Palsy.   In  section  5,  the  recursive  partitioning  classification 
rule  given  in  Friedman  (1977)  is  used  to  construct  a  classification 
scheme  based  on  the  groupings  obtained  by  the  clustering  method.   This 
classification  scheme  is  useful  in  examining  the  effects  of  corrective 
measures  such  as  surgery  or  assistive  devices  on  the  gait  of  individual 
patients.   The  results  given  in  Section  5  indicate  that,  in  general, 
operation  is  effective  while  the  significance  of  assistive  devices  are 
hard  to  detect. 
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2.   DESCRIPTION  OF  DATA 
The  sample  used  for  this  study  is  consisted  of  62  patients  at 
Boston's  Children's  Hospital,  ranging  in  age  from  three  to  twenty- 
four,  all  of  whom  were  dignosed  as  having  some  of  the  traits  associ- 
ated with  Cerebral  Palsy,   Using  high-speed  movie  cameras,  these 
patients  were  filmed  walking.   These  films  were  then  digitized  and 
translated  into  measurements  associated  with  leg  motion.   Measurements 
used  for  the  analysis  were  taken  over  the  course  of  one  gait  cycle. 
A  gait  cycle  begins  when  the  heel  or  toe  of  one  foot  strikes  the 
floor,  and  ends  when  that  same  heel  or  toe  strikes  the  floor  next. 
Measurements  associated  with  the  speed  of  gait,  the  length  of  steps, 
the  time  spent  on  one  leg  or  both  legs,  and  the  angular  motion  of 
the  hip,  knee,  and  ankle  were  available  for  all  patients.   Those  used 
included: 

*cycle  time — how  many  seconds  to  complete  a  gait  cycle. 
*cadence — the  number  of  steps  per  minute,  based  on  the  cycle 

time  observed. 
*swing  time — time  spent  with  only  one  foot  on  the  floor; 

measured  as  a  percentage  of  cycle  time  for  both  right 
and  left  sides. 
*double  limb  stance  time — time  spent  with  both  feet  on  the 

floor;  measured  as  a  percentage  of  cycle 
time  for  both  right  foot  forward  and  left 
foot  forward. 
*step  length — measured  for  both  steps. 

*stride  width — distance  in  inches  between  heels  in  the  frontal  plane. 
*cycle  velocity — feet  per  minute,  based  on  cycle  time  and  step  length. 


-4- 


The  following  measurements  of  angular  motion  were  made  for  the  left  and 
right  sides  of  each  patient: 

*hip  motion,  x-direction — movement  in  the  sagittal  plane,  taken  from  the 
side  view.   The  measurement  was  arrived  at  by  finding  the 
angle  between  the  line  segment  for  the  sacral  stick  (a 
stick  extending  straight  out  from  the  patient's  back)  and 
the  line  segment  from  the  hip  center  to  the  knee  center.   If 
the  angle  measured  90  degrees,  hip  motion  in  this  direction 
was  defined  to  be  zero  degrees.   See  (a)  in  Figure  A. 

*hip  motion,  y-direction — movement  in  the  frontal  plane.   It  is  measured 
from  the  angle  formed  by  the  line  segment  between  the  marker 
spots  on  the  pelvic  belt  and  the  line  segment  between  the  hip 
and  knee  centers.   If  this  angle  is  90  degrees,  then  motion  in 
this  direction  is  defined  to  be  zero  degrees.   See  (b)  in 
Figure  A. 

*hip  motion,  z-direction — motion  in  the  transverse  plane  which  requires 

readings  to  be  taken  from  the  front  and  side  views  at  corres- 
ponding times  of  the  cycle.   From  each  view,  a  line  segment 
was  found  which  intersected  the  front  of  the  knee,  and  was 
perpendicular  to  the  line  between  the  hip  and  ankle  centers. 
The  angle  of  rotation  is  measured  to  be  zero  if  the  hip, 
patella  and  ankle  fall  in  a  straight  line  as  seen  from  the 
front.   Movement  of  the  patella  away  from  the  hip  is  defined 
as  rotation.   See  (c)  in  Figure  A. 
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*knee  motion,  x-direct ion — measured  from  the  angle  formed  by  the  line 
segments  from  the  hip  to  knee  centers  and  the  knee  to 
ankle  centers.   Zero  degrees  motion  occurs  when  the  two 
line  segments  fall  in  a  straight  line.   See  (d)  in  Figure  A. 
*ankle  motion,  x-direction — measured  from  the  angle  formed  by  the  line 

segment  from  the  knee  center  to  ankle  center  and  the  line  of 
the  floor.   See  (e)  in  Figure  A. 
*ankle  motion,  z-direction — measures  foot  rotation.   The  measurement  points 
are  taken  to  be  the  center  of  the  ankle  and  a  marker  placed 
between  the  second  and  third  toes.   When  these  two  points 
line  up  from  the  front,  the  angle  of  rotation  is  zero  degrees. 
See  (f)  in  Figure  A. 
(Remarks:   It  should  be  pointed  out  that  measurements  of  foot  rotation 
are  much  more  subject  to  measurement  error  than  the  other  angular  measure- 
ments because  of  the  small  distance  between  measurement  points.   And 
foot  rotation  measurements  for  the  patients'  left  side  were  not  included 
in  this  analysis  due  to  large  amounts  of  missing  data) . 

Each  of  the  angles  described  above  was  measured  in  each  frame  of 
film  throughout  a  gait  cycle.   (The  cameras  were  run  at  a  speed  of 
fifty  frames  per  second.)   The  resulting  measurements  form  a  time 
series   x  ,  which  can  be  approximately  represented  by  a  thirteen  terms 
Fourier  series  x(t)  =  A^  +  A^  cos(t)  +  B  sin(t)  +  ...  +  A^cos(6t)  +  B^ 
sin(6t)   for  each  frame  t,   using  the  least  squares  method.   For  each 
patient,  the  set  of  coefficients  so  obtained  for  a  given  angle  is  used 
to  describe  the  corresponding  angular  motion. 
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Because  the  thirteen  terms  representation  of  each  angle  resulted 
in  an  excessive  number  of  variables  available  for  analysis,  a  pro- 
cedure was  required  which  would  significantly  reduce  the  number  of 
coefficients  used  for  each  angle.   The  informal  rule  that  was  used 
for  selecting  coefficients  consisted  of  two  steps.   In  Step  One,  for 
any  given  patient  and  for  any  angular  motion,  those  coefficients 
which  were  over  twice  as  large  in  magnitude  as  the  root  mean  square 
(RMS)  were  defined  as  "significant"  coefficients;  that  is,  coeffici- 
ent C  is  significant  only  if   |c(/RMS  ^  2,  where 

N 


RMS  = 


E   (x  -  5  )V(N  -  13) 
i=l 


X   is  the  fitted  value  of  the  series  at  frame  t,   and  N  is  the 
number  of  frames  used  to  film  the  gait  cycle.   This  definition 
was  then  applied  to  a  random  subset  of  twelve  patients  in  Step 
Two,  and  a  coefficient  was  included  in  the  analysis  if  two  or 
more  patients  had  "significant"  values  for  that  coefficient.   As 
a  result  of  this  procedure,  37  terms  are  used  to  represent  the 
eleven  different  angular  motions  instead  of  the  orginal  143 
possible.   Although  it  lacks  the  formal  significance  of  t-tests, 
this  method  of  data  reduction  is  not  an  unreasonable  one  as  the 
high  noise  levels  of  the  angular  measurements  render  them  not 
to  be  too  useful  for  fine  discrimination  between  patients  anyway. 
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3.   CLUSTER  ANALYSIS  OF  GAIT  MEASUREMENTS 
One  purpose  of  this  study  is  to  group  the  sampled  Cerebral 
Palsy  patients  on  the  basis  of  the  various  gait  measurements, 
with  the  hope  that  the  resulting  clusters  would  identify  with 
different  severeness  levels  of  abnormal  gait.   Since  the  gait 
measurements  of  a  patient  can  be  perceived  as  a  point  in  the 
46-dimensional  Euclidean  space,  and  clusters  of  patients  can 
be  thought  of  as  high-density  point-clouds  in  this  measurements- 
space  separated  from  other  such  point-clouds  by  regions  of  low 
point-density,  the  high-density  clustering  model  given  in 
Hartigan  (1975,  p.  205)  and  Wong  (1981)  is  appropriate  for 
this  clustering  problem.   Under  this  clustering  model,  the 
kth  nearest  neighbour  clustering  procedure  given  in  Wong  (1981) 
is  known  to  be  asymptotically  set-consistent  for  high-density 
clusters;  it  has  also  been  shown  to  be  useful  for  identifying 
high-density  clusters  in  small  samples.   Hence  this  clustering 
technique  is  used  here  to  obtain  a  hierarchical  clustering  of 
the  Cerebral  Palsy  patients. 

The  clustering  tree  obtained  when  the  kth  nearest  neighbour 
clustering  procedure  (with  k  =  2)  is  applied  to  the  gait  measure- 
ments of  the  62  patients  is  shown  in  Figure  B.   This  sample 
hierarchical  clustering  corresponds  to  the  tree  of  high-density 
clusters  defined  on  the  kth  nearest  neighbour  estimate  of  the 
underlying  measurements-space  density.   Five  modal  clusters  can 
be  identified  in  the  sample  clustering  tree,  but  two  of  these 
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clusters  (modes  3  and  5,  as  shown  in  Figure  B)  may  correspond  to 
local  bumps  present  in  the  estimated  density  function  when  k  =  2   is 
used;  indeed,  for  larger  values  of  k(k  =  3,  4,  and  5),  only  two 
modal  clusters  (modes  1  and  4)  can  be  identified  in  the  result- 
ing hierarchical  clusterings.  However,  since  it  is  of  clinical 
interest  to  explore  and  find  any  homogeneous  subgroup  of 
patients  whose  gait  characteristics  are  different  from  the  majority 
of  Cerebral  Palsy  patients,  the  clustering  solution  obtained 
for  k  =  2  will  be  further  examined  in  the  following  section, 
in  spite  of  the  possibly  artificial  nature  of  some  of  the 
identified  subpopulations .   The  labels  of  the  patients  in  each 
modal  cluster  shown  in  Figure  B  are  listed  below: 
Cluster  1: 


Cluster  2 


Cluster  3 


Cluster  4 


Cluster  5 


S473.R7  S561.R1 

S506.R16  S499.R2 

S455.R3  S497.R6 

S439.R3  S413.R4 

S427.R3  S428.R4 

S470.R3  S507.R4 


S405.R10 

S593.R6 

S452.R4 

S454.R1 

S429.R11 

S542.R8 


S592.R3 
S450.R4 
S400.R2 
S398.R6 
S432.R4 
S415.R4 


S595.R9   S410.R3 
S557.R8 

S568.R1 
S434.R3 
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4.   CLUSTER  PROFILES  AND  CHARACTERISTICS 

Summary  statistics  of  the  five  identified  clusters  were  computed  for 
each  of  the  46  gait  measurements,  and  they  are  given  in  Table  1,   More- 
over, for  each  measurement,  a  one-way  analysis  of  variance  was  performed 
across  clusters  and  the  resulting  F-ratio  is  also  shown  in  Table  1. 
Since  the  groupings  are  specif iced  by  the  clustering  procedure,  the  com- 
puted F-ratio  will  not  conform  to  the  usual  F-distribution  but  will  be 
larger  just  by  chance,  so  they  cannot  be  used  to  perform  the  usual  signi- 
ficance testing.   However,  these  pseudo-F  ratios  serve  the  purpose  of 
identifying  individual  variables  that  are  very  different  across  the 
five  clusters;  for  example,  cycle  velocity  (F  =  115.9)  is  seen  to  be  dras- 
tically different  among  the  clusters. 


Insert  Table  1  here 


To  illustrate  and  compare  the  characteristics  of  the  various  clus- 
ters, graphical  profiles  of  the  mean  vector  of  the  strength  measurements 
(Figure  C)  the  hip  movements  (Figure  D) ,  the  knee  movements  (Figure  E) , 
and  the  ankle  movements  (Figure  F)  for  each  cluster  are  plotted.   In  all 
of  these  figures,  the  measurements  are  identified  on  the  horizontal  axis. 
In  Figure  C,  the  vertical  axis  shows  the  relative  magnitude  of  the  mean 
for  each  cluster  to  the  largest  cluster  mean  for  that  variable  (or  measure- 
ment).  Cluster  3  is  distinguished  from  the  rest  by  having  longer  cycle  time. 
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shorter  swing  time  and  longer  double-limb  stance  time,  shorter  step  lengths, 
and  lower  cycle  velocity.   Patients  in  Cluster  1  appear  to  be  physically 
stronger  than  those  in  Cluster  3,  but  they  are  consistently  worse  than 
those  in  the  remaining  clusters.   Clusters  2,  4,  and  5  exhibit  similar 
characteristics,  having  relatively  short  cycle  time,  fast  cadence,  long 
swing  time,  and  high  cycle  velocity;  however.  Cluster  2  can  be  seen  to 
be  the  strongest  group  physically.  Hence,  Clusters  3,  1,  4,  5  and  2  are 
in  the  order  of  increasing  gait  strength. 

The  absolute  relative  magnitudes  of  the  means  of  the  hip  movement  co- 
efficients are  shown  in  Figure  D.   The  corresponding  plots  for  knee  and 
ankle  movements  are  given  in  Figure  E  and  Figure  F  respectively.   It 
should  be  pointed  out  that  the  absolute  value  of  each  Fourier  coefficient 
reflects  its  importance  in  determining  the  overall  shape  of  the  original 
waveform;  hence  the  larger  the  absolute  value  of  the  coefficient,  the  more 
effect  it  has  in  determining  the  shape  of  the  waveform.   For  the  hip  and 
ankle  movements  (Figures  D  and  F) ,  it  can  be  seen  that  the  absolute  rela- 
tive magnitudes  of  the  Fourier  coefficients  for  Cluster  3  tend  to  be  very 
small,  indicating  less  angular  motion  in  the  hip  and  akle  than  the  other 
clusters.   On  the  other  hand,  Clusters  2,  4,  and  5  consistently  have  large 
absolute  relative  magnitudes  in  these  movements,  indicating  a  higher  degree 
of  muscular  control  over  these  leg  motions,  and  hence  a  less  pathological 
gait  (see  also  the  relative  physical  strength  shown  by  these  clusters  in 
Figure  C) ;  it  can  also  be  noted  that  Cluster  4  distinguishes  itself  from 
the  other  two  clusters  by  showing  relatively  limited  right  ankle  movements 
(see  Figure  F) .   However,  no  systematic  pattern  across  clusters  can  be  detected 
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for  the  knee  coefficients  shown  in  Figure  E.   And  the  above  observations 
point  to  the  significance  of  hip  and  knee  movements  in  indicating  the 
physical  strength  of  a  patient's  leg  motion,  while  showing  that  the  cor- 
relation between  the  knee  movements  and  the  gait  strength  of  a  patient 
is  very  low. 

The  variation  of  medical  history  of  the  patients  across  clusters 
are  also  examined,  and  are  given  in  Table  2.   Most  of  the  patients  in  the 
two  physically  weakest  groups  (Clusters  1  and  3)  were  diagnosed  as  having 
quadriplegia  (i.e.  all  four  limbs  are  affected),  while  none  of  the 
patients  in  the  physically  strongest  group  (Cluster  2)  had  the  same  diag- 
nosis.  On  the  other  hand,  a  patient  diagnosed  as  having  diaplegia  (i.e.  only 
two  limbs  are  affected)  can  be  in  any  one  of  the  five  identified  clusters; 
hence,  the  present  diagnosis  system  is  not  suitable  for  indicating  the  physical 
status  of  a  Cerebral  Palsy  patient.   Moreover,  it  can  be  seen  in  Table  2 
that  the  patients  in  Clusters  1  and  3  were  more  likely  to  have  been  operated 
on  (69%)  than  the  other  three  clusters  (38%).   It  is  apparent  that  physicians 
are  presently  more  inclined  to  recommending  surgery  to  physically  weaker 
patients  than  to  the  stronger  patients. 


Insert  Table  2  here 
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5.   DISCUSSION  AND  FURTHER  ANALYSIS 
One  purpose  of  this  study  is  to  identify  the  different  severity 
stages,  if  any,  of  Cerebral  Palsy  by  examining  the  gait  measurements  of 
a  sample  of  Cerebral  Palsy  patients.   In  Sections  3  and  4,  it  is  shown 
that  several  severity  stages  of  Cerebral  Palsy  could  be  identified  using 
the  kth  nearest  neighbour  clustering  procedure;  gait  cycle  velocity  and 
other  strength  characteristics  play  a  predominant  role  in  distinguishing 
the  identified  stages,  although  the  magnitudes  of  hip  and  ankle  movements 
are  also  useful  in  discriminating  the  various  stages.   However,  the  validity 
of  the  identified  clusters  is  still  questionable,  and  additional  samples 
are  needed  for  further  analysis  and  cross-validation  purposes.   Alternatively, 
the  discriminating  power  of  the  clustering  procedure  used  can  be  determined 
by  examining  its  effectiveness  in  distinguishing  samples  from  known  distinct 
populations  with  different  gait  patterns,  but  the  relevant  gait  data  have 
yet  to  be  measured  or  collected. 

Subsequent  to  cluster  analysis  which  is  useful  for  groupings-identification, 
it  is  important  to  develop  a  classification  scheme  so  that  patients  can  be 
assigned  to  one  of  the  identified  clusters.   Such  a  classification  scheme 
might  provide  a  useful  typology  for  classifying  Cerebral  Palsy  patients;  for 
example,  it  might  be  found  that  the  effectiveness  of  corrective  measures, 
such  as  surgery  or  assistive  devices  like  crutches,  varies  across  the  clus- 
ters (or  equivalently,  across  different  severity  stages) . 

For  illustration  purposes,  the  recursive  partitioning  decision  rule  given 
in  Friedman  (1977)  is  used  to  construct  a  classification  scheme  based  on  the 
groupings  obtained  by  the  clustering  method.   The  resulting  classification 
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tree  is  given  in  Table  3.   By  virtue  of  his  or  her  gait  measurements, 
a  patient  can  then  be  classified  into  any  one  of  the  five  identified  clusters 
by  descending  along  the  appropriate  branches  of  the  classification  tree. 
Among  the  gait  measurements  considered  in  this  study,  the  most  important  one 
in  discriminating  between  the  identified  clusters  is  seen  to  be  gait  cycle 
velocity,  which  is  used  for  discrimination  along  most  branches  of  the  class- 
ification tree.   However,  the  angular  motion  of  the  right  knee  is  important 
in  distinguishing  Cluster  4  from  Clusters  2  and  5. 

Although  the  error  rate  found  in  this  classification  scheme  is  only  1 
out  of  29,  it  is  a  deflated  figure  and  is  not  useful  for  evaluating  the 
effectiveness  of  the  classification  procedure  because  the  original  group- 
ings are  based  on  the  results  of  a  clustering  method.   But  this  classifi- 
cation tree  is  useful  in  assigning  a  new  gait  observation  to  one  of  the 
identified  clusters.   In  three  different  instances,  measurements  taken  on 
a  patient  before  and  after  surgery  resulted  in  different  assignments  to  the 
identified  groups  according  to  the  above  classification  scheme;  in  one 
case,  the  patient  improved  from  Cluster  A  to  Cluster  2  after  knee  surgery, 
and  in  the  remaining  two  cases,  the  patients  moved  up  from  Cluster  3  to 
Cluster  1.   But  similar  improvements  were  not  detected  for  patients  who 
were  given  assistive  devices  such  as  braces. 


Insert  Table  3  here 
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Table  1.   SUMMARY  STATISTICS  for  CLUSTERS 


Clust 

er  1 

Clust 

er  2 

Cluster  3 

Cluster  4 

Cluster  5 

Measurements 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

F-Ratlo 

Cycle  Time- 
Cadence 

1.04 

0.19 

0.81 

0.07 

2.05 

0.31 

0.86 

0.08 

0.90 

0.11 

51.3 

119.13 

22.73 

150.00 

13.61 

59.63 

9.66 

139.71 

11.45 

134.79 

15.06 

25.4 

L.  Swlnj;  7, 

35.14 

3.42 

40.67 

2.66 

20.87 

11.05 

37.15 

3.64 

37.03 

3.54 

11.6 

K.  Doul)li;  % 

15.92 

4.19 

10.19 

1.75 

27.31 

8.00 

12.22 

3.73 

12.02 

3.03 

13.6 

L.  Double  % 

13.75 

3.81 

11.47 

1.05 

36.30 

13.20 

10.35 

1.08 

13.70 

3.77 

18.8 

R.  Step  ].i:nf) 
L.    Step  Leng 

th 

1 .  18 

0.23 

1.60 

0.11 

0.67 

0.29 

1.54 

0.14 

1.48 

0.19 

19.1 

th 

1 .  29 

0.37 

1.65 

0.10 

0.60 

0.47 

1.66 

0.14 

1.57 

0.22 

12.2 

Stride  Width 

0.41 

0.11 

0.47 

0.10 

0.55 

0.27 

0.44 

0.19 

0.43 

0.10 

0.8 

Cycle  VelocI 
K.  Hip  X 

ty 

143.05 

13.94 

243.39 

18.25 

39.12 

21.81 

223.98 

30.92 

203.15 

12.22 

115.9 

Ao 

5.83 

5.01 

0.34 

6.04 

-0.09 

9.13 

-20.00 

9.76 

6.56 

3.18 

17.2 

Ai 

-3.54 

2.83 

-2.55 

2.96 

-0.09 

2.28 

11.00 

5.67 

-6.54 

3.90 

21.4 

Bi 

-3.32 

3.86 

-6.21 

4.45 

-1.73 

3.44 

-4,98 

5.91 

-3.54 

4.96 

0.7 

R.  Hip  Y 

Ao 

-6.95 

14.04 

-28.03 

15.40 

-15.07 

10.21 

17.78 

14.06 

-10.83 

18.70 

5.8 

Ai 

1.33 

5.06 

-3.17 

1.34 

0.49 

4.53 

3.80 

9.50 

0.48 

2.39 

1.1 

R.  Hip  Z 

Ao 

43.93 

10.45 

30.67 

8.79 

40.86 

15.24 

-34.38 

7.39 

34.38 

4.73 

63.2 

Ai 

13.32 

6.25 

6.99 

4.49 

1.13 

8.69 

-10.02 

9.06 

12.09 

5.40 

12.9 

Az 

0.11 

3.63 

1.03 

1.86 

0.49 

3.54 

1.83 

3.80 

-1.97 

5.51 

0.9 

Bi 

9.57 

4.93 

19.04 

4.74 

5.89 

5.37 

-13.65 

6.43 

16.91 

10.22 

18.2 

B2 

-7.04 

5.09 

-10.29 

8.71 

-3.01 

3.05 

-7.26 

4.95 

3.99 

10.88 

4.4 

1,.  Hip  X 

Ao 

6.33 

4.93 

3.24 

2.86 

2.11 

4.48 

-10.81 

11.55 

2.67 

3.45 

8.0 

Ai 

5.87 

3.23 

6.89 

4.66 

0.80 

2.38 

-5.15 

4.98 

5.38 

3.15 

11.3 

Bi 

1.94 

3.21 

3.30 

5.91 

2.23 

1.21 

0.44 

10.28 

2.34 

6.88 

0.1 

L.  Hip  Y 

Ao 

-12.11 

14.45 

9.62 

5.38 

-6.81 

10.39 

11.04 

24.82 

-18.41 

15.90 

4.5 

Ai 

-2.67 

5.76 

2.57 

3.98 

-1.38 

4.21 

-3.47 

5.77 

3.16 

2,96 

3.1 

L.  Hip  Z 

Ao 

37.17 

3.89 

37.46 

6.80 

44.06 

13.97 

-28.04 

10.33 

31.43 

5.79 

72.5 

Ai 

-9.50 

6.17 

-5.03 

5.45 

-1.09 

6.57 

6.45 

3.07 

-10.83 

6.28 

9.3 

Az 

2.63 

2.62 

3.41 

1.82 

0.88 

3.29 

1.17 

4.77 

2.36 

3.56 

0.5 

Bi 

-12.20 

5.50 

-21.35 

4.30 

-3.52 

5.83 

15.23 

5.34 

-17.38 

7.27 

31.5 

B2 

-4.66 

3.83 

-7.05 

4.77 

0.89 

3.64 

-4.76 

5.57 

4.89 

12.69 

3.0 

R.  Knee  X 

Ao 

-8.94 

7.96 

-15.04 

8.29 

-6.13 

11.50 

20.18 

8.00 

-11.96 

9.97 

12.3 

Ai 

0.71 

3.45 

-0.91 

2.20 

-1.09 

1.75 

-1.45 

2.90 

-0.05 

3.32 

0.7 

A2 

0.87 

1.70 

0.30 

3.75 

0.88 

1.38 

1.06 

1.99 

-0.54 

2.72 

0.7 

Bi 

1.66 

2.89 

-4.54 

8.41 

-0.16 

2.30 

2.63 

5.87 

-1.23 

6.15 

1.6 

B2 

-0.21 

2.42 

1.40 

3.48 

0.80 

1.06 

-1.07 

2.36 

-0.11 

2.94 

0.7 

L.  Knee  X 

Ao 

-16.43 

9.81 

-1.47 

2.43 

-7.76 

7.12 

15.56 

16.63 

-18.55 

8.45 

12.4 

Bi 

0.13 

4.46 

-4.10 

3.70 

-0.30 

2.78 

-1.81 

11.72 

4.50 

4.84 

2.1 

B2 

0.29 

2.01 

-0.10 

2.13 

-0.81 

1.57 

-0.04 

2.59 

0.72 

4.08 

0.3 

R.  Ankle  X 

Ao 

8.44 

13.18 

12.63 

11.21 

-2.36 

20.22 

-0.47 

8.31 

26.21 

18 .  69 

4.2 

Ai 

0.46 

3.38 

3.43 

3.87 

-0.70 

5.21 

-0.91 

5.88 

4.56 

4.16 

2.4 

Bi 

-0.36 

5.30 

-0.82 

5.15 

0.39 

1.75 

0.80 

4.50 

-1.39 

5.27 

0.2 

R.  Ankle  /, 

Ao 

32.41 

9.10 

33.82 

5.28 

40.70 

7.39 

-29.29 

7.31 

28.62 

10.76 

54.7 

Ai 

8.66 

4.96 

17.10 

5.20 

2.33 

5.63 

-11.41 

7.62 

8.85 

5.47 

17.7 

Bi 

-16.60 

5.90 

-20.63 

2.62 

-6.51 

8.07 

16.14 

3.93 

-21.67 

5.82 

41.0 

L.  Ankle  X 

Ao 

19.70 

10.30 

10.64 

10.59 

7.17 

23.40 

-1.63 

2.91 

38.67 

8.37 

10.5 
3.4 

Ai 

-0.54 

3.91 

-5.96 

3.80 

-1.59 

4.27 

-3.07 

3.88 

-6.01 

3.57 

'  *  i 

Bi 

-0.22 

2.73 

-0.25 

3.70 

-1.16 

2.95 

-1.91 

5.61 

0.58 

2.16 

0.6 
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Table  3.   Recursive  Partitioning  Classification  Tree 


N 

-  29 

Cluster 

1 

2 

3 

4. 

5 

f 

11 

4 

5 

5 

4 

Cycle  Velocity  s  174.79 


N  <=  16 

Cluster 

1 

2 

3 

4 

5 

// 

11 

0 

5 

0 

0 

Cycle 

Velocity  5  95.89 


N  -  5 

Cluster 

1 

2  3 

4 

5 

f 

0 

0  5 

0 

0 

Cycle 

Velocity  >  95.89 


N 

-  11 

Clusteri 

1 

2 

3 

4 

5 

# 

11 

0 

0 

0 

0 

Cycle  Velocity  >  174.79 


N  =  13 

Cluster 

\ 

2 

3 

4 

5 

« 

0 

4 

0 

5 

4 

\/ 


Klghf  Angle  Z 
Ao  £  17.66 


N  =  5 

Cluster 

1 

2 

3 

4 

5 

// 

0 

0 

0 

5 

0 

ghC  Angle  Z 
Ao  >  17.66 


H 

N  -  8 

Cluster 

1 

2 

3 

4 

5 

tf 

0 

4 

0 

0 

4 

Cycle 
Velocity 
S  217.23 


Cycle 
Velocity 
>  217.23 
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-FIGURE  A       ANGULAR  MEASUREMENTS  OF  LEG  MOTION 
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